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Abstract. Many heavy fermion materials are known to crossover from the Kondo 
lattice regime to the mixed-valent regime or vice-versa as a function of pressure or 
doping. We study this crossover theoretically by employing the periodic Anderson 
model within the framework of the dynamical mean field theory. Changes occurring 
in the dynamics and transport across this crossover are highlighted. As the valence 
is decreased (increased) relative to the Kondo lattice regime, the Kondo resonance 
broadens significantly, while the lower (upper) Hubbard band moves closer to the Fermi 
level. The resistivity develops a two peak structure in the mixed valent regime: a low 
temperature coherence peak and a high temperature 'Hubbard band' peak. These two 
peaks merge yielding a broad shallow maximum upon decreasing the valence further. 
The optical conductivity, likewise exhibits an unusual absorption feature (shoulder) 
in the deep mid-infrared region, which grows in intensity with decreasing valence. 
The involvement of the Hubbard bands in dc transport, and of the effective /-level 
in the optical conductivity arc shown to be responsible for the anomalous transport 
properties. A two-band hybridization-gap model, which neglects incoherent effects 
due to many-body scattering, commonly employed to understand the optical response 
in these materials is shown to be inadequate, especially in the mixed-valent regime. 
Comparison of theory with experiment carried out for (a) dc resistivities of CcRhlns , 
Ce2Ni3Si5, CeFeCea and YbIr2Si2; (b) pressure dependent resistivity of YbInAu2 
and CeCue; and (c) optical conductivity measurements in YbIr2Si2 yields excellent 
agreement. 
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1. Introduction 

Rare earth lanthanides and actinides [1, 2] exhibit a wide range of behaviour such 
as heavy fermions (HFs), mixed valence (MV), proximity of superconductivity and 
magnetism, quantum critical points etc. Such behaviour arises through an interplay 
of a variety of factors such as hybridization between conduction bands and deep /- 
levels, orbital degeneracy, crystal field effects, long range spin interactions and most 
importantly local Coulomb repulsion [3-7]. 

In this work, we present a detailed and systematic theoretical investigation of a 
regime that borders on heavy fermions at one end and on the mixed valence regime 
at the other. Experimentally, such a crossover from HFs to MV or vice-versa has 
been observed to happen through pressure or doping [8-11]. The effects of such a 
crossover have been investigated for several materials [12-21]. For example, the first 
known heavy fermion superconductor CeCu2Si2 [9], when doped with Yttrium (a non- 
magnetic homologue of Ce) shows increasing mixed valent character. The resistivity 
exhibits a two peak structure in the temperature range 2K to 300K. Varying doping 
concentration results in gradual coalescing of the two peak structure into a single broad 
peak. For the heavy fermion compound YbCu4Ag [10], the ambient pressure resistivity 
as a function of temperature shows a broad peak which is characteristic of mixed valent 
compounds. By applying pressure, the broad peak sharpens and the peak position also 
shifts to lower values of T, thus indicating a crossover to the Kondo lattice regime. 
The coefficient of the low temperature Fermi-liquid resistivity also increases sharply 
with pressure, implying a decrease in the coherence scale. CeBeis [11] also shows a 
pressure-induced crossover from Kondo lattice regime to mixed- valent regime, as seen 
in the changes in the coherence peak in the resistivity. The trend in this material is 
opposite to that seen in the previous example, YbCu4Ag. [10] 

Most previous theoretical attempts to describe the effects of pressure on heavy 
fermion materials have employed the single-impurity Anderson model [22-24]. An 
illustrative and important work in this context is that of Chandran et al [22]. They 
use a phenomenological model comprising a competition between elastic energy cost 
and valence fluctuations induced magnetic energy gain. The magnetic energy for an 
Anderson lattice model was computed by using the free energy of an impurity Anderson 
model and ignoring lattice coherence effects. The free energy itself was arrived at 
through a slave-boson mean field approximation [25] which is a static approximation 
and is thus unable to treat dynamical effects of valence fluctuations. Phenomenological 
expressions were used to model volume dependence of the parameters, and the pressure 
was obtained by using P{T,V) = —dF/dV. The authors were able to describe 
continuous and discontinuous valence transitions in a single framework. However, since 
the impurity Anderson model was used, lattice coherence effects were ignored. Transport 
quantities were not calculated. Recent years have seen the use of lattice models to 
describe the concentrated Kondo systems. The minimalist model that accounts for 
a large part of heavy fermion and mixed-valent behaviour is the periodic Anderson 
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model (PAM) which represents a lattice of localized /-orbitals with a Hubbard repulsion 
U hybridizing locally with a wide non-interacting conduction band. The PAM has 
of course been explored extensively using a wide range of methods and techniques 
such as dynamical mean field theory (DMFT) [26,27], extended DMFT [28], cluster 
expansions [29], finite size simulations [30] etc. Within the framework of dynamical 
mean field theory (DMFT) [26,27], the self-energy becomes purely local or momentum 
independent, which simplifies the problem significantly, while retaining the competition 
between intinerancy and localization in a non-trivial way. The momentum independence 
of the self-energy implies that we can view lattice problems as locally self-consistent 
impurity problems, hence the solution of PAM, for instance reduces to that of a single- 
impurity Anderson model with a self-consistent hybridization. The various impurity 
solvers that have been adapted to solve the effective impurity problem are numerical 
renormalization group [31], quantum Monte Carlo [32], exact diagonalization [33], 
perturbation theory methods such as iterated perturbation theory [34], local moment 
approaches [35], and slave-particle approaches [36] etc. 

The Kondo lattice regime within the PAM has been investigated heavily by 
us [35,37,38] and the other groups [31,39,40] within dynamical mean field theory 
and cluster extensions ignoring the d— f repulsion effects. We have used local moment 
approach(LMA) within the DMFT framework previously to understand the PAM in 
the Kondo lattice (KL)limit [37]. The LMA is a non-perturbative diagrammatic theory 
based approach developed by Logan and co-workers [41]. This approach has been 
benchmarked extensively against methods such as Bethe Ansatz [41] and numerical 
renormalization group [42]. Excellent quantitative agreement has been found, thus 
providing justification for its use as an impurity solver within DMFT. Further advantages 
of using LMA are, that real frequency quantities are obtained directly with reasonable 
computation expense at all temperature and interaction strengths. Our focus in 
the LMA+DMFT approach to PAM was on universality and scaling in dynamics 
and transport. A single low energy scale was found to characterize the spectra and 
transport in the Kondo lattice regime [35]. The mixed- valent regime was not our focus, 
nevertheless, our studies indicated the absence of universality and scaling, although 
adiabatic continuity to the non-interacting limit was seen. The dc and optical transport 
properties were compared to several Kondo insulators and heavy fermion metals and 
excellent quantitative agreement was found [38,43]. A few mixed- valent materials 
such as YbAla and the skutterudite compound CeOs4Sbi2 were also considered [44] 
and again good agreement between theory and experiment was found. In a recent 
work, valence transition in Ytterbium and Europium intermetallics was studied by 
Zlatic and Freerics [45] employing a multi-component Falicov-Kimball (FK) model 
within dynamical mean field theory. The authors argue that a complete description 
would entail a solution of the periodic Anderson model combined with the FK model. 
However, since this is challenging, they choose to solve just the FK model, albeit a multi- 
component one. Transport quantities like dc and optical conductivity, thermopower and 
magnetoresistance etc were calculated and qualitatively compared to the experiment. 
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Using an equation of motion decoupling approximation, Bennard and Coqblin [46] 
have investigated the PAM and explored the variation of the /-valence with various 
parameters of the model. They use the results of this study to understand pressure 
dependent valence changes. Miyake and co-workers [47] have used a variety of methods 
including the slave-boson approximation, the density matrix renormalization group 
etc to investigate the one dimensional extended PAM (EPAM), which includes the 
d — f repulsion effects represented by Ucf- The EPAM has been investigated within 
dynamical mean field theory by Sugibayashi and Hirashima [48]using quantum Monte- 
Carlo methods. Their focus has been to understand the interplay between Uc/, valence 
fluctuations and superconductivity. In a recent work Ylvisaker [49] et. al combined 
local density approximation with dynamical mean field theory to understand the valence 
fluctuation and the valence transition in the Yb metal. For the impurity solver, they 
have used Hirsh-Fye quantum Monte-Carlo and continuous time quantum Monte-Carlo. 
They reproduce the experimentally observed valence transition, and conclude that the 
Yb metal is a fluctuating valence material rather than an intermediate-valent one. As 
mentioned before and as illustrated through the above mentioned studies, there has 
been substantial work on the valence transition and its effect on the spectral quantities 
within the PAM, however the effects on transport quantities due to valence fluctuations 
and the crossover regime between the Kondo lattice and the mixed-valent regime have 
received scant attention. 

In this work, we focus on such a KL-MV crossover using the LMA-I-DMFT approach 
to the PAM. We assume that effects of pressure/doping would be to change the model 
parameters, and hence a scan of the parameter space within the PAM framework should 
be able to provide insight into the crossover regime. We highlight the changes occuring 
in the dynamics and transport properties as a result of this crossover. We find several 
new results such as a two peak resistivity and anomalous absorption features in the 
optical conductivity in certain parameter regimes. We provide theoretical explanations 
for these anomalies, and show that such behaviour does indeed exist in real materials and 
may be explained quantitatively using the present approach. The paper is structured 
as follows: We present the model and formalism in the next section. The results and 
discussions are in sections 3 and 4 respectively. A comprehensive range of experimental 
measurements is shown to be described by our theoetical results in section 5. We 
conclude in section 6. 

2. Model and Formalism 

The periodic Anderson model is one of the simplest models representing a paradigm 
for understanding the physics of heavy-fermion compounds. In standard notation, the 
Hamiltonian for PAM is given by 
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The first two terms describe the c-orbital energy (ec) and the kinetic energy of 
the conduction band arising from nearest neighbour hopping t, which is scaled as 
t oc it* is the unit of energy) in the large dimension limit where the coordination 



number Z^. — t- oo. We choose the hyper cubic lattice for our calculation, for which the 
non-interacting density of states given by Po(f) = exp (— e^/t^)/t*-y/7r is an unbounded 
Gaussian [26]. The hybridization between the c- and /-electrons is represented by the 
third term {V) and is responsible for making otherwise localised /-electrons itinerant. 
The last term is /-orbital site energy (e/) and the on-site Coulomb repulsion U for two 
/-electrons of opposite spin. 

In the non- interacting {U = 0) case, the two orbital energies ec and e/ completely 
determine the nature of the ground state, i.e. whether the system would be gapped or 
gapless at the Fermi level. Equivalently, the occupation numbers nf = '^ZaUlafia) ^^^1 
Uc = (4(7^10-) s-lso be used to characterize the ground state. It is straightforward 
to see that the system is a band insulator for Uf + Uc = 2, while it is metallic for 
Uf + He 7^ 2 [4]. In the absence of hybridization, i.e. ^ = 0, the / and c sites decouple 
from each other. Since the /-part is site-diagonal in this limit, the /-Green functions 
are trivial to obtain, e.g. through the equation of motion method [50]. The tight- 
binding part along with the on-site c-orbital energy may be solved through Feenberg 
renormalized perturbation theory (FRPT) [50], and the non-interacting Green's function 
may be obtained as gQ^uj) = {uj~^ — Cc — S{uj))~^. The Feenberg self-energy S{u) is in 
fact a lattice-specific functional of g^, For example, 5'[(?q] = jt^gQ for the Bethe lattice 
in the limit of infinite dimensions [50]. 

Within DMFT, the /-selfenergy, is local [26,27] in real space, hence the FRPT 
may be used for V ^ also to obtain the / and c-Green's functions. These are then 
given by [35,43] 

UJ — 6/ — ^{UJ 



CO 



uj^ — e^ — S{uj) 



— — S{uj) 



(2) 



(3) 



a;+ — e/ — 

where ^(0;) is the momentum- independent /-self-energy representing interaction effects 
and 5'(ci;) is the same functional of G'^ as it is of g^ in the non-interacting limit. This 
Feenberg self-energy may be obtained through the relations given below (please see 
references [35,43] for more details): 

G\u) = H[^] = r = , , , (4) 

where -^[7] denotes the Hilbert transform of 7 with respect to the non-interacting 
density of states po{e). While 7 = in the ^ = limit, for finite V, it is given 

by 7(0;) = 00+ -ec- V^[oo+ -€f- S(a;)]-i. 
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2.1. Local moment approach: theory and practice 

As mentioned in the introduction, within DMFT, the lattice model is mapped onto an 
impurity model with a self-consistently determined hybridization. And the solution of 
the resulting impurity problem is the most challenging step in the solution of the lattice 
model. We choose the local moment approach within DMFT to solve the PAM. The 
LMA has been benchmarked extensively in case of the SIAM with other exact methods 
such as NRG and Bethe Ansatz. While in case of the PAM, excellent agreement with 
experiments for a wide variety of materials ranging from Kondo insulators to heavy 
fermion metals and for a range of experiments such as dc and optical transport as well 
as magnetotransport justifies its use to understand the crossover from Kondo lattice 
regime to the empty orbital regime via the mixed valence regime. The theoretical 
formalism and practical implementation are described in references [35,37,43], and the 
reader is referred to these works for details. 



2.2. Transport: linear response 



Within DMFT, vertex corrections are absent, hence the single-particle Green's functions 
are sufficient within Kubo formalism to obtain transport quantities such as dc resistivity 
and optical conductivity. The expressions derived previously [37] can be brought into a 
simpler form (derived in the appendix) given below: 



a{uj) = 0Re 



°° iripiu') — Upioj + u') 
duj 



GHuj + io') GHuj') - GHuj 



7(0; + w') — 7*(u;') 7(0; + w') — 7(0;') 

where cxo = Ane'^t'^a'^n/h for lattice constant a, electronic charge e, and electron density 
n and 7(0;) = — — V'^[uj~^ ~ ^ ^('^)]~^- The DC conductivity is obtained by 
considering the w — limit (see appendix for more details) of the above equation as 



27r2 



Re 



du 



drip 
du 



TiD\u) 
Im7(u;) 



2(l-7(a;)G») 



(6) 



where D'^{u) = —lmG'^{u)/'n' is the spectral function of the retarded Green's function 



G'^(u) and upiu) 



1) ^ is the Fermi function. 



(5) 



3. Results 



In this section, we show results for spectral functions, i.e the density of states (dos), 
at both zero and finite temperature, optical conductivity again for T > 0, and DC 
resistivity. We work with a reasonably large interaction strength {U ^ 4.5t*) and 
a hybridization strength of = O-Gt^. The crossover from Kondo lattice regime 
to the mixed- valent regime is investigated by varying ec and e/ such that the total 
occupancy Utot = Uc + Uf is fixed, amd n/ is decreased. This is done because pressure 
experiments would be expected to keep the total occupancy fixed. The practical details 
of implementation are given in references [35,37]. We begin with the low energy scale 
and the spectra. 
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3.1. Coherence scale 

The low energy coherence scale has been identified in previous studies to be ul = 
ZV"^ /t^, where Z = {1 — (9S/(9a;|^=o)^^ is the quasiparticle weight (inverse effective 
mass). This scale is exponentially small in the strong coupling Kondo lattice regime 
{uf — 7- 1). Spectral functions, optical conductivities and resistivity were shown to be 
universal functions of {T /ujl.uj /ul) in our previous studies [35,37,43]. Figure 1 shows 
the variation of the lattice coherence scale ujl with the /-orbital occupancy Uf for 
fixed Utot = 1-25 (filled circles) and 1.1 (squares). A crossover from the Kondo lattice 
regime {uf — j- 1) through the mixed-valence regime (0.3 ^ Uf < 0.8) to the empty 
orbital regime(n/ <^ 1) manifests in a rapid increase in ul- A lowering of the /-orbital 




0.4 0.6 0.8 1 



Figure 1. A sharp decrease of the lattice coherence scale, wl is seen, with n/ 
increasing from the empty orbital regime (n/ 1) through the mixed- valence regime 
(0.3 < n/ < 0.8) to the Kondo lattice regime (?i/ — > 1). 

occupation implies greater local charge fiuctuations, implying an effective decrease of 
correlation effects. This in turn implies a decrease in effective mass or an increase in Z 
and hence an increase in ul with decreasing n/. 

3.2. Density of states: Zero temperature 

The density of states or the spectral functions are considered next. Changes in the low 
frequency Kondo resonance and in the high frequency Hubbard bands are identified. 

We show the T = spectral functions {D^{uj; T = 0) = -ImG^(w,T = 0)/7r u = 
f,c) as a function of the absolute frequency u/t^, in the left panel of figure 2. In the 
/-dos (solid line), a usual three peak structure with high energy Hubbard bands and a 
narrow Kondo resonance at the Fermi level is seen for higher nf{> 0.7). For lower nj, 
a new spectral feature (marked by an arrow) appears above the Fermi level, the weight 
of which grows with decreasing Uf. The lower Hubbard band (LHB) moves closer to 
the Fermi level, while the upper Hubbard band (UHB) shifts to higher energies with 
decreasing uj. The middle panel in figure 2 shows the Kondo resonance in greater 
detail as a function of the scaled frequency u/ul. Close to the Kondo lattice regime, 
a pseudogap is seen straddling the resonance, which gets progressively filled up and for 
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the lower values of rij, there is no trace of a pseudogap. Such a transfer of spectral 
weight is very non-trivial and is a very significant feature since it manifests clearly in a 
second peak in optical conductivity, as will be discussed in sections 3.4 and 4.2. 
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Figure 2. (colour online) Left panel: The zero temperature spectral functions, 
D"{u}\T = 0), I' = c(dashed), /(solid) are shown as a function of 'absolute' frequency 
w/t* for various /-occupancies (nj, indicated in the boxes and in the legends) and a 
fixed ritot = 1-25. The middle and the right panels show the same spectra as the left 
panel, but as a function of the scaled frequency a; /w^ highlighting the changes in the 
Kondo resonance and the pseudogap. 



The conduction band density of states, D^{{jJ]T = 0), shown in the left panel (red 
dashed lines) of figure 2, has an overall Gaussian envelope, with spectral weight carved 
out at the effective /-level, (seen clearly in the right panel of figure 2 and further 
discussed in section 4.1) in the form of a pseudogap and small Hubbard bands flanking 
the envelope. In contrast to the /-dos (left panel of figure 2), the Hubbard bands here 
possess a very small fraction of the total spectral weight. The lower Hubbard band is 
distinct from the envelope at higher nf(> 0.6), and merges into the envelope at lower 
r2/(< 0.5). The pseudogap fills up with decreasing rif and for the lowest rif values 
shown, there is indeed no trace of a pseudogap (see inset). The position of the Kondo 
resonance, the Hubbard bands and a few other aspects are discussed in greater detail 
in subsection 4.1. 



3. 3. Density of states: Temperature dependence 

In figure 3, we show the evolution of the /-dos as a function of temperature for nj = 0.8 
and Uf = 0.3 in the left and right panels respectively. The insets in the panels show the 
low frequency, low temperature {{u/ujl,T/ujl) ~ part more clearly. 

The left inset shows that the pseudogap proximal to the Fermi level fills up, while 
a thermal broadening of the Kondo resonance is clearly visible in the /-dos. At any 
given temperature T, it is expected that changes in the spectral function will occur in 
the frequency range w < T. And indeed, this is seen in the KL regime as discussed 
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Figure 3. The /-electron spectral function at temperatures T = T /lol = 0, 2, 5, 10, 20 
(indicated in legends) and for all scales in the main panel. The insets show the low 
frequency part as a function of scaled frequency for lower temperatures. The left panel 
is for Uf = 0.80 and the right panel is for Uf = 0.30 (with fixed ntot = 1.25). 

below and in reference [37]. The main left panel shows spectra for three different 
temperatures, which are hard to distinguish. This implies that there has hardly been 
any spectral weight transfer on a scale of t*, even at the highest temperature shown 
(T = 10, equivalent to ~ O.lSt*). 

The physics in the mixed- valent regime is, naturally, different. The coherence scale 
is not exponentially small (see figure 1). In this regime, spectral weight transfer occurs 
in an energy interval that is far greater than the thermal energy scale. This may be seen 
in the right panel of figure 3, where at a temperature of 5cJi(~ 0.77t,,), the UHB which 
is at ~ 6t* and hence about 15 times the temperature gets strongly affected. In fact, 
the integrated spectral weight in the UHB is found to decrease with T. Such differences 
in terms of spectral weight transfer between the KL and MV regime have been noted 
previously [37]. The new spectral feature above the Fermi level (marked by an arrow 
in the left panel of figure 2) that emerges distinctly in the MV regime gains spectral 
weight and grows with increasing T. The reason for the growth of this feature is not 
completely clear; however, since the total spectral weight is conserved, the melting of 
Kondo resonance and the loss of spectral weight in the UHB with increasing T must be 
compensated by gain in spectral weight elsewhere, and this could be one of the reasons 
for the growth of this feature with increasing T. 

Angle-resolved photoemission experiments should be able to easily identify such a 
feature. The conduction band density of states exhibits similar temperature dependence 
as the /-dos and is hence not shown. 
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3.4- T = optical conductivity 

As mentioned before, the zero temperature Green's functions are sufficient within DMFT 
to compute the T = optical conductivity, a{uj] T = 0) (equations 5) [26,27]. In figure 4, 
we show a{uj\T = 0) as a function of frequency {oj/t^) for various values of and a 
fixed Utot = 1.1. For rij — 1 , a familiar structure of a single mid-infrared peak (MIR) 
arising due to interband transitions [51], marked with a thin dashed arrow in figure 4, is 
obtained. As we decrease n/, the MIR peak shifts to higher frequencies. A Drude peak 
is obtained at the lower frequencies for all n^, but it cannot be shown here since it is a 
Dirac delta function for T = 0. 




co/t^ co/t^ 

Figure 4. The zero temperature optical conductivity (t(u;;T = 0) as a function of 
frequency (w/t*) for various n/, and fixed ntot = 1-1 (left panel) and 1.25 (right 
panel). The MIR peak is marked with a thin dashed arrow, while the DMIR feature 
is marked with a thick solid arrow. 

The lower nj values hold a surprise. A deep mid- infrared (DMIR) absorption feature 
or a shoulder, marked with a thick solid arrow in figure 4, emerges for lower Uf values. 
This is distinct from the above mentioned MIR peak. Such a two peak structure has 
been reported earlier for the classic mixed valent compound CePds [52] and recently 
for YbIr2Si2 [12]. For the former compound, spin-orbit splitting was argued to be 
responsible for the two peaks, but the latter has not been understood quantitatively. 
Recently it has been pointed out through band structure studies that the DMIR feature 
could be due to transitions to the effective /-level [53] . Here we confirm that this feature 
could indeed arise from absorption into the effective /-level at = Z{ef + S(0)), at 
which frequency, a pseudogap appears in the spectral functions (figure 2). Thus the 
scenario proposed here could be relevant for both CePds and YbIr2Si2. A detailed 
discussion of the effective /-level, the pseudogap and transfer of spectral weight into the 
pseudogap region is given in section 4.2. 
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Here we must add that a two-peak structure is indeed obtained for higher rifoj, such 
as rttot = 1-25 shown in the right panel of figure 4 and hence is a generic feature in 
the crossover regime. Furthermore, with decreasing n/, the DMIR feature merges into 
the MIR peak, without being able to develop into a full distinct peak for lower nj, and 
hence gets harder to distinguish for higher ritot- 

We now move on to finite temperature transport. 



3.5. DC Resistivity as a function of temperature 

In figure 5, we show the dc resistivity as a function of temperature T/t*. The left panel 
shows the resistivity for 0.65 < nj < 0.95, while the right panel is for 0.30 < Uf < 0.60. 
The theoretically computed p(T) all have zero residual values, but for visual clarity, 
they have been appropriately vertically offset (mentioned in the caption). There is a 
clear shift of the coherence peak (low temperature peak) to higher temperatures since 
its position correlates with the low energy scale [37], and ul itself increases sharply with 
decreasing Uf. For nf < 0.8, a second maximum begins to form at higher temperature. 
This second peak become easily distinguishable in the mixed- valence regime and merges 
into the coherence peak to yield a broad maximum close to the empty orbital regime (see 
right panel, Uf < 0.45). Here we have shown that the crossover from the Kondo lattice 
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Figure 5. DC resistivity for various nj values (mentioned next to the curve) as a 
function of temperature T/i* for a fixed n^ot = 1.25. The theoretically computed 
p(r) all have zero residual values, but for visual clarity, they have been appropriately 
vertically offset. For example - in the left panel, we add a temperature independent 
po that is zero for ri/ = 0.65 and increases by unity for every higher rij (po — 1, 2, 3... 
for n/ = 0.7,0.75,0.8,...). 



regime to the mixed-valent regime involves the coalescing of the coherence peak and a 
high temperature peak into a single broad peak. This high temperature peak is not 
present in the KL regime or the empty orbital regime, and appears in the mixed-valent 
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regime only due to the proximity of the Hubbard band to the Fermi level (discussed in 
detail in section 4.3). Although we have displayed data in figure 5 for ritot = 1.25, the 
two peak structure is a generic feature of the mixed- valent regime for other ntot as well. 
However, the degree to which the two peaks can be distinguished and identified varies 
in different parameter regimes. 

Pressure/doping-dependent resistivity measurements in CeCu2Si2 [8,9,20] show 
similar behaviour. For low pressure, a two peak resistivity is observed. As a function 
of increasing pressure, the two peaks merge and a single broad peak is seen. In the 
experimental paper, the effect of pressure was ascribed to a crossover from Kondo lattice 
to intermediate valence, however the second peak was argued to arise from higher lying 
crystal-field multiplets. Here, although we do concur with the crossover phenomenon, 
we present an alternative scenario for a two peak resistivity structure. The scenario is 
that with increasing pressure, the valence decreases, bringing the Hubbard band in the 
proximity of the Fermi level, thus leading to a second peak arising due to the Hubbard 
band states contribution. We emphasise here that a two-peak structure is obtained 
in our theory despite the model Hamiltonian (equation 1) not containing crystal-field 
effects. 

We now move on to finite temperature optical transport. 
3.6. T > Optical transport 




(O/t* co/t* 



Figure 6. Optical conductivity (7(0;; T) for Uf — 0.75, ric = 0.61 and rif ^ 0.40, ric = 
0.71 in the left and right panels respectively (ec = 0.5t, and U ^ 4.7t,) as a function 
of w/t, in the main panels and as a function of w/ojl in the insets respectively for 
various temperature values (given as multiples oi ujl in the legends). The insets show 
cr(a;; T) for low temperatures to show the transfer of spectral weight more clearly. 

We now show the temperature dependence of the optical conductivity in figure 6 
for two different /-occupancies, 0.75 and 0.4 (other parameters are mentioned in the 
figure caption) in the left and right panels respectively. The a{u;T) for the higher Uf 
case has a usual prominent mid-infrared (MIR) peak, with a slight shoulder at lower 
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frequencies. The inset shows the low temperature evolution of (T(a;), while the main 
panel shows T < lOa;^. With increasing temperature, a transfer of spectral weight over 
all energy scales in seen. The lower n/(0.4) optical conductivity shown in the right 
panel of figure 6 is qualitatively different from that of the left panel. The shoulder like 
structure seen in the Kondo lattice regime (left panel) develops into a full peak in the 
mixed- valent regime, while the MIR peak diminishes in intensity and gets broader. The 
analysis at T = shows (see section 4) that the lower frequency peak occurs at the 
effective /-level e*j. The inset shows that the DMIR peak or the shoulder appears only 
below a temperature T = Q.QSujl even though it is positioned at a frequency u equal 
to the coherence scale ul (since here ej + S(0) ~ 1). This precise behaviour is seen in 
recent optical conductivity measurements in YbIr2Si2 [12] (see section 5). 

4. Discussion 

4-1 ■ Spectral features: Kondo resonance and Hubbard bands 

In this subsection, we discuss the positions and weights of the various features seen in 
the spectra in section 3.2. At high energies, the prominent spectral features are the 
Hubbard bands. In figure 2, it is seen that the Hubbard bands shift to the right with 
decreasing nj. This is expected from the atomic limit, since in that limit, the Hubbard 
bands occur at and ef + U, and ej moves closer to the Fermi level when nj decreases. 
However, we find from our calculations, that in the displayed spectra, the Hubbard bands 
are not at the positions predicted by the atomic limit, as must be expected, due to a 
combination of shifts in the levels arising through hybridization and self-energy effects. 
We see that even for Uf ~ 0.8, the LHB and the UHB are not symmetrically placed 
about the Fermi level as would have been the case in the strong coupling KL regime. In 
fact, the LHB is at w < — t*, while the UHB is at cu > U. The LHB is predicted very well 
by the unrestricted Hartree Fock (UHF) solution (not discussed here) which considers 
the static part of the self energy and hybridization effects. For predicting the UHB, we 
can neglect hybridization effects, however, we need to retain the static part and most 
importantly the real part of the self energy. The latter may be approximated by an 
expression obtained through high frequency moment expansion [54], and thus the UHB 
turns out to be at ~ {ef + VepT^) /2, where ef = ef + U{n + fi) /2, tl^ = U^n{2 - fi) 
and have been defined in section 4.2 of reference [35] as the fictitious occupation 
number and the moment obtained from the spin-dependent host/medium propagators 
that are used to construct the LMA self-energies. 

In the other extreme, i.e close to the Fermi level, a low frequency form of the 
Green's functions using a Fermi liquid self-energy S(u;) = S(0) + ~ 1) + C>{u^) may 
be derived that enables us to understand most of the features of the Kondo resonance. 
The resulting spectral functions in the neighbourhood of the Fermi level are given by 
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where e*j = Zef = Z(e/- + S(0)) is the effective /-level These equations show that the 
low energy {u < ul) spectral features are precisely those of a non-interacting {U = 0) 
PAM with renormalized parameters [V"^ — )■ ZV^, ef — )■ Z{ef + S(0)) and (ec — )■ €c) [35]. 
Thus, this is the renormalized non- interacting limit (RNIL). A scaling collapse of 
the numerically obtained spectra with the analytical expressions above would be a 
demonstration of adiabatic continuity of the interacting system to the non-interacting 
limit. Indeed, we do see such a collapse (not shown), implying adiabatic continuity, 
although the range of frequencies over which such a collapse occurs decreases from 
w < in the KL regime, to a; ^ in the empty orbital regime. 

As C(j — )■ in equations 7, the renormalized non-interacting limit spectral functions 
also vanish, thus the pseudogap seen in the spectra is positioned at the effective /- 
level, e*j = Z (e/ -|- S(0)). Transfer of spectral weight into the pseudogap happens 
with decreasing as seen in figure 2 which implies that the frequency interval in 
the neighbourhood of e*j- is gaining spectral weight. This spectral weight transfer is 
completely missed by theories that are equivalent to the renormalized non-interacting 
limit. This includes approaches such as slave-boson theories or two-band models of 
heavy fermion systems, which ignore the imaginary part of self-energy. The variation of 
e*jr with Uf may be easily predicted using the Luttinger's theorem which states [35] 

_(n^ + = y po(e) de + ^(-e/) . (8) 

For fixed ritot = nf + Uc, the upper limit of integration on the right side of equation 8 
is also fixed. Crossing over from the KL to MV regime would require decreasing rib- 
and increasing n^, which in turn would require decreasing Cc. Thus to keep — + V^/ej 
fixed, the e/ = e/ + S(0) must increase, and hence the effective /-level, = Zef, also 
increases as Uf decreases. This is indeed seen in the figure 2 because the position of the 
pseudogap does indeed shift to higher frequencies as ri/ is decreased. 



4-.2. Inadequacy of the renormalized non-interacting limit 



In [37], it was shown that the two band model or the RNIL predicts a square root 
singularity at the minimum direct gap (~ 2\/ZV), and hence the MIR peak is generally 
attributed to the direct gap. A simple Fermi-liquid analysis T,{u) ^ ReS(0)+a;(l — 1/Z) 
of the poles of the k-dependent conduction electron Green's function 



G"^(c^;ek) = 
yields a two-band model 



(ec + ek - e}) ± J (e^ + ek - e)) + AZV^ 



(9) 



(10) 
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The minimum direct gap is given by min(aj+(ek) — a;_(ek)). The square root singularity 
at the direct gap appears in the two band model because the imaginary part of the 
self energy is neglected, and the resulting spectral functions are Dirac delta functions. 
Including incoherent effects due to electron-electron scattering results in broadening of 
the MIR peak and cutting off the square root singularity. 

At the lowest values, a clear two-peak structure is visible in the optical 
conductivity displayed in the figure 4 (especially for rif = 0.45 in the left panel). The 
low frequency peak is in fact at the effective /-level, as argued below, while the high 
frequency peak is the usual MIR peak. Naively, finding an absorption peak at is 
counter-intuitive, because the renormalized non-interacting limit (equations 7) shows 
that a pseudogap exists at e^, implying that there is no density of states at that energy. 
So how can absorption into a gap happen? The answer is of course that the RNIL, 
which is equivalent to a slave-boson mean-field theory, which in turn is equivalent to 
a two-band model [51], are not totally correct in their predictions. These approaches 
neglect the imaginary part of the self-energy (scattering rate). So even though the RNIL 
predicts a gap, there is in fact no gap at when self-energy effects are included (see the 
pseudogap feature in the density of states in figure 2, and the discussion in section 4.1). 
And in fact, it may be shown rigorously, without recourse to LMA that absorption to 
e*j will happen provided the imaginary part of self energy is non-zero at that energy. 
To see this, consider the roots of the real part of the denominator of the conduction 
electron Green function, G'^(a;,ek) (equation 9), given by 

Re [cj — ec — ek — V'^ [u — ef — S(cj))~"'"] = 

oj-ec-V^ ^ = Ck (11) 

If the imaginary part of the self energy is completely neglected and a first order Taylor 
expansion is carried out for the real part of the self energy, we get back equation 10. 
However, retaining the imaginary part, however small it might be, results in an equation 
that is at least cubic in order. And for ek = Cc — Z{ef + S(0)), one of the roots is just 
e} = Z(e; + S(0)). 

We support the arguments above with LMA results below. The dispersion, aj(ek) 
is computed for rif = 0.95(nc = 0.58), 0.7(?ic = 0.63) and Uf = 0.4(nc = 0.71) with 
= 0.5, U ~ 4.7t=i, through equation 11 and shown in the top, middle and bottom panels 
respectively in figure 7. The two-bands obtained at the renormalized non-interacting 
level (eq. 10) are also superimposed in red. It is clear from the figure that the agreement 
of the two-band model with the full dispersion gets progressively worse as n/ decreases. 
For rif = 0.4, the middle band is clearly visible in the full dispersion, while being 
completely absent in the two-band picture. This third band, as argued above is centred 
at a; = e*f. It is the excitations from the 'band' below the Fermi level to e^, that 
appears as an additional 'anomalous' absorption peak in the deep mid-infrared region 
(as a shoulder). This anomalous peak is seen to become prominent in the mixed- valent 
regime and is very small or invisible in the Kondo lattice regime. We refer the reader 
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Figure 7. The dispersion a;(ek) computed from equation 11 (black) and from the 
RNIL, equation lO(red) for = 0.5i*, U ^ 4.7t, and nj = 0.95, ric = 0.58 (top panel), 
rif — 0.70, ric — 0.63 (middle panel ) and rif — 0.40, ric — 0.71 (bottom panel). 

to figure 4, where an extra absorption feature is easily discernible for lower Uf values, 
while for higher Uf values, just a weak shoulder is observed. As we will see later, 
such a two-peak optical conductivity has indeed been observed in recent experiments 
(sections). 

4-3. Two-peak resistivity: why? 

We have shown previously that the coherence peak (low temperature peak) in the 
resistivity, which would be a minimum in the conductivity [37], occurs at a temperature 
comparable to ul, the low temperature scale. A second peak (high temperature peak) 
has indeed been observed in experiments and has been attributed to crystal field split 
levels(for example in Ref [8]). Notwithstanding the foregoing possibility, we find an 
alternative and far simpler explanation for the existence of the second peak even within 
the single band PAM. The second peak occurs at a temperature that is roughly half 
of the lower Hubbard band energy scale. It might at first seem surprising to note that 
the Hubbard band is contributing to transport. The Hubbard bands are usually at an 
energy scale of f//2 in the strong coupling Kondo lattice regime which being of the order 
of a few eV, remain untouched until room temperature. Nevertheless, for mixed-valent 
systems, even if the U is large, either the LHB or the UHB moves close enough to the 
Fermi level so as to be affected at room temperature scales. 

To show the contribution of the Hubbard band states to the conductivity, we 
consider expression 6 again. On the right hand side, the integration is carried out over 
all frequencies. However, if we introduce a upper and a lower cutoff in the integration 
limits, we can isolate the contribution of individual spectral features within those limits. 
In figure 8, we show the calculated dc resistivity for rij = 0.6, = 0.65 (same as that in 
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9- 




T/(0, 



Figure 8. Figure shows that the second peak in the DC resistivity seen most clearly in 
the right panel of figure 5 for n/ = 0.6 disappears progressively by excluding the lower 
hubbard band and other spectral features. The solid line is for = 0.6, ric = 0.66 
and is the full resistivity. The dashed and the dotted lines are obtained by excluding 
the lower Hubbard band below the Fermi level and with decreasing cutoff above the 
Fermi level (see text for discussion). 

figure 5), with three different cutoffs: (i)the sohd hne being no cutoff (full resistivity), 
(ii) the dashed hne having hmits such that the lower Hubbard band is excluded but 
the new 'non-Hubbard band' feature above the Fermi level (marked with an arrow in 
figure 2) is included and (iii) the dotted line having limits such that both the LHB 
and the new feature is excluded. It is seen that excluding the spectral weight of the 
Hubbard band and the new feature enhances the resistivity systematically. This shows 
that the Hubbard band contribution to the dc conductivity is substantial in mixed-valent 
compounds. The two-peak behaviour of resistivity, is not specific to the parameter 
regime for which the dc resistivity has been displayed {U/V"^ ~ 8,ntot = 1.25). In other 
parameter regimes, such as ritot = 1-1, our calculations demonstrate (not shown here) 
that the dc resistivity crosses over from a single sharp coherence peak to a broad peak 
lineshape with decreasing valence, going via a two-peak structure but with very faint or 
hard to distinguish features. 

5. Comparison to Experiments 

5.1. DC resistivity: ambient pressure 

In figure 9, we have superposed theoretically computed dc resistivity (solid lines) 
with the experiments (circles) for CeRhlns [55], Ce2Ni3Si5 [56], CeFeGes [57] and 
YbIr2Si2 [12]. To find the appropriate theoretical parameters for each material, we 
had to adopt a trial and error approach. However, a surprising similarity in the 
dc resistivity of CeRhlus and Ce2Ni3Si5 reduced our effort substantially. We found 
that we could take the experimentally measured dc resistivity for the two materials 
and scale them onto each other simply by rescaling the x and y axes. CeRhlns 
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has anisotropic resistivity, nevertheless the pa and pc may also be scaled onto each 
other, thus showing that qualitatively, they are also similar. We find that these 
two materials have an /-occupancy ~ 0.7, with the rest of the parameters being 
Uc ~ 0.55, f//V^ ~ 8,ec = 0.5t*. For CeFeGes, we found that the mixed-valence 
resistivities do not fit the data well. Rather, the best fit was found using parameters 
{jif ~ l,r2c ~ 0.77, f//y^ ~ 5, ec = 0.3) that signify an intermediate correlation 
with /-occupancy being nearly unity. The optical transport data for YbIr2Si2 (see 
below) shows a two-peak structure very similar to that seen in figure 4 for nj = 0.4 
and ritot = l-l- Thus we take the same resistivity and superimpose that onto the 
experimentally measured one, and we see very good agreement. 

50 
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Figure 10. Comparison of pressure dependent resistivity measurements (magnetic 
contribution only) for YbInAu2 with theory. The inset shows the reported [58] 
experimental data. 
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5.2. DC resistivity: Pressure dependence 

Hydrostatic pressure dependence measurements of the dc resistivity of YbInAu2 and 
its non-magnetic homologue LuInAu2 have been carried out by Fuse et.al [58]. The 
experimental magnetic resistivity along with the residual resistivity {pmag + Po) is shown 
in the inset of figure 10. The main panel shows a comparison of theory (solid lines, 
Utot = 1-1,U/V'^ ~ 8) with experiment (symbols). The agreement is seen to be excellent. 
The Tmax is seen to decrease with pressure, and the material appears to be progressing 
towards a Yb^"*" state with increasing pressure, as conjectured in the experimental 
paper. Nevertheless, the valence remains in the mixed- valent regime (< 0.4), even 
with pressures upto 20KBar. 

We now move on to CeCug, which has gathered a lot of attention in the past decade 
as a material that can be tuned to a quantum critical point with Au doping. Pressure 
dependent resistivity measurements on CeCug were carried out in 1985 by Thompson 
and Fisk [59]. They found that with increasing pressure, the material crosses over 
from Kondo lattice like to mixed- valence like regime. In the left panel of figure 11, we 
show the experimental graph, while in the right panel, the theory (same data as that for 
figure 5) is shown. The experiment shows that with increasing pressure, the temperature 
dependence follows the P = curve to some temperature and then deviates. In the inset 
of the left panel, the experimental data for < P < 17.4kbar is shown to collapse when 
plotted as R/Rmax vs. T/Tmax for T < iTmax- The theoretical curves in the right panel 
correspond very closely to those seen in the experiment. The agreement between theory 
and experiment shows that indeed with increasing pressure, the occupancy does change 
from Uf ^1 to njr^ 0.75, thus implying that the valence of Ce changes from ~ 3+ to 
~ 2.75+. The theory curves decrease more rapidly (than experiment) with increasing 
temperature beyond the coherence peak, and the reason for this is that the phonon 
contribution is not subtracted in the experiment. The inset in the theory panel has five 
different valencies plotted together, namely, nj = 0.90,0.85,0.80,0.78 and 0.75. Except 
for Uf = 0.75, the rest of the data is seen to collapse onto a single curve, as seen in the 
experiment. The uj = 0.75 curve does collapse upto a certain temperature and then 
deviates from the universal curve. 

5.3. Optical transport: YbIr2Si2 

The recently discovered heavy fermion system YbIr2Si2 [60] has a crystal structure 
similar to the well studied YbRh2Si2. The latter exhibits a field-tuned quantum critical 
point (QCP), while the former has a pressure-tuned first order phase transition to a 
ferromagnetic phase. The experimentally measured optical conductivity [12] is shown in 
the top panel of figure 12. A clear two peak structure is evident, and a large scale spectral 
weight transfer occurs as temperature is increased from 0.4K to 300K. An important 
characteristic of the temperature dependence of a{u; T) is that, as temperature is 
increased, the lower frequency peak (shoulder) merges into the continuum at T ~ 60K, 
at which temperature the higher frequency peak remains untouched. 
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Figure 11. Left panel: Experimentally measured pressure dependent resistivity 
measurements for CeCug [59]. Right panel: Theoretically computed p{T) (same data 
as left panel of figure 5) plotted as p/ Pmax vs- T/Tmax- The respective valencies are 
mentioned next to the curves. 



We note that the theoretically computed optical conductivity shown in figure 6 
does appear to resemble the experimentally measured a{u]T). In the experiment, 
the shoulder peak appears only for T < 30K, which, according to theory, should be 
roughly 0.08 times the peak frequency for the chosen parameters(see the discussion 
for figure 6). Thus, theoretically, we can predict that the shoulder should appear 
at ~ 30.0/0.08/^ ~ 32meV. The shoulder peak position as predicted by theory 
indeed agrees very well with the experimentally observed shoulder position (~ 30meV). 
The bottom panel of the same figure shows the optical conductivity computed for 
U/t^ ~ 5, ~ O.Gtl and ec = O.St^,, which yields rif = 0.4 and = 0.7, thus classifying 
YbIr2Si2 as a mixed- valent material. These parameters were chosen using the results 
shown in figure 4, and for consistency, we note that the dc resistivity in the experiment 
has a broad and shallow peak (see bottom right panel of figure 9) which indicates that 
YbIr2Si2 belongs to the mixed- valent regime. The DC conductivity obtained through a 
low frequency limit in the theory is higher than the corresponding experimental values. 
This is natural, since the theory neglects electron-phonon scattering, which if included 
would reduce the theoretical conductivities. Apart from this disagreement, an excellent 
agreement between theory and experiment is seen in quantitative terms for both the 
lineshape and the temperature dependence. A phenomenological analysis [12] of the 
experimental optical conductivity using Drude and extended Drude formalism has been 
used to infer non-Fermi liquid nature of the quasiparticles in this material. However, 
if the quantitative agreement between theory and experiment is any indication, this 
material is a perfect Fermi liquid. The self energy has the correct Fermi liquid form 
for the parameters used to compute the optical conductivity in the bottom panel of 
figure 12. In our approach, the Luttinger's theorem is used as a constraint, which 
also supports the inference of Fermi liquid behaviour. The low temperature resistivity 
also exhibits a clean behaviour, both in the experiment [13] and theory (shown to 
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Figure 12. top panel: Experimentally measured optical conductivity for YbIr2Si2. [12] 
Bottom panel: Theoretically computed cr(a;;r) for U /t^, ^ 5,V'^ ^ 0.6^^ and 
Cc = 0.5i*, which yields rif = 0.4 and ric = 0.7. Excellent quantitative agreement 
is observed between theory and experiment for the line shape and the temperature 
dependence. 



agree well in the inset of the bottom right panel in figure 9), which again supports a 
Fermi liquid ground state. Additionally, we conclude that this material is in the mixed- 
valent regime with an effective /-occupancy of n/ = 0.4, which would correspond to 
an effective valence Yb^'^^. (Naturally, the Uf given here must be interpreted as hole 
occupancy for Yb compounds.) This effective valence predicted by our theory disagrees 
with the inference from susceptibility measurements, which yields a Yb valence close to 
3+. Further theoretical and experimental investigations are needed to resolve this. We 
can also infer that the effective /-level given by e*j = Z{ef + S(0)) in YbIr2Si2 is located 
at roughly 25-30 meV. The dc resistivity calculated from theory for the parameters 
mentioned in the figure 12 does agree qualitatively with experiment in terms of a broad 
peak and the lineshape(see figure 9), however, the scale inferred from fitting to the 
experimental data does not agree with that obtained from optical conductivity. There 
appears to be an inconsistency in theory, in terms of an ambiguity in the parameters that 
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can fit the experiment. The implication is, we believe, minor. It should be possible by a 
better scanning of the parameter space to find a 'right' set of parameters that can fit the 
optical and dc conductivity consistently with the same scale. We believe that finding this 
right set of parameters will surely make the theory more convincing, but probably not by 
a great measure, since almost all aspects of the experiment are quantitatively captured 
by the presently used set of parameters. The new set of parameters would probably 
change the numbers (e.g for the [/, V"^ values), nevertheless, the major conclusions such 
as that of YbIr2Si2 being a Fermi liquid and a mixed- valent system should remain 
unscathed. 

6. Conclusions 

We have carried out an extensive study of the changes in spectral and transport 
properties of the periodic Anderson model as the /-occupation number is varied from 
unity to nearly zero. The local moment approach within dynamical mean field theory 
has been employed for this study. We have used the results of our study to understand 
the crossover from mixed-valent regime to the Kondo lattice regime (or vice-versa) 
observed to occur in many rare earth intermetallics. We find several unusual features 
such as a two peak resistivity, anomalous absorption feature in optical conductivity etc 
in the crossover regime. We show that the proximity of the Hubbard band to the Fermi 
level is responsible for the former while the latter happens due to optical transitions 
into the effective /-level. The two-band (hybridization-gap) model generally applied 
to understand the physics of these materials is shown to be inadequate in this regime. 
Qualitative agreement with pressure-dependent dc resistivity measurements in CeCu2Si2 
is found while quantitative agreement with dc transport measurements in CeRhlns, 
Ce2Ni3Si5, CeFeGes and YbIr2Si2 is obtained. We also obtain excellent agreement with 
pressure dependent resistivity measurements in YbInAu2 and the prototypical quantum 
critical system CeCug. In agreement with previous observations, we find that increasing 
pressure pushes Cerium materials to mixed- valence while the Yb materials cross over to 
the Kondo lattice regime. Further, a remarkable agreement with optical conductivity 
experiments in YbIr2Si2 is obtained. We also infer from the agreement that YbIr2Si2 
belongs to the mixed-valent regime. Further investigations including d — f correlations 
in the periodic Anderson model are underway to understand valence transitions. 
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Appendix 

The expression for dynamical conductivity (hypercubic lattice) in [37] is 

/•oo /"OO 

T) = — dui [nF{ui) - npicoi + co)] / de po{e)D'{e; UJi)D''{e; u + uji) (12) 

where D'^{e;u) = —lmG'^{e;u})/'iT;G'^{e;uj) = (7(0;) — e)""*^ and 7(0;) is defined below 
equation 4. The identity lmzilmz2 = Re{zlz2 — ZiZ2)/2 for zi, Z2 & C may be used to 
get 



27r2 \^ 7(wi) - 7*(a; + Wi) 7(^1) - 7(0; + Wi) 

The Hilbert transform integral over e in equation 12 may be carried out to give the 
expression 5. The — )■ limit of the equation 5 yields the dc conductivity. The 
second term involves using the L'Hopital's rule which in turn requires the knowledge 
of dG'^{uj)/d'y, which we derive below. The c-Green's function is given by a Hilbert 
transform over the Gaussian dos. The result is expressed in terms of the complementary 
error function as [26] 

G'{uj) = H[-f] = -is0Fexp(-72)erfc(-is7) (13) 

where s = sgn{lm'y) = +1 for the retarded functions considered here. The derivative of 
this with respect to 7 is straightforward and is given by 

= -2jG%u) - ^sv^exp(-7')(-^^)^^^l„, = 2 (1 - ^G%u)) (14) 
Using this yields equation 6. 
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